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Abstract 

The low frequency magnetic field spectrum in the primordial plasma is of 
particular interest as a possible origin of magnetic fields in the universe (e.g., 
Tajima et al. 1992 and Cable and Tajima 1992). We derive the magnetic field 
spectrum in the primordial plasma, in particular, at the epoch of primordial 
nucleosynthesis. The pioneering study of Cable and Tajima (1992) of the 
electromagnetic fluctuations, based on the Fluctuation-Dissipation Theorem, 
is extended. Our model describes both the thermal and collisional effects in 
a plasma. It is based on a kinetic description with the BGK collision term. 
It is shown that the zero-frequency peak found by Cable and Tajima (1992) 
decreases. At high frequencies, the blackbody spectrum is obtained naturally 
without the necessity of the link procedure used by them. At low frequencies 
(u) < 4w pe , where u) pe is the electron plasma frequency) it is shown that the 
magnetic field spectrum has more energy than the blackbody spectrum in 
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vacuum. 

PACS numbers: 52.25.Dg, 52.25.Gj, 98.80.Ft 
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I. INTRODUCTION 



Although plasma is the main constituent of the primordial universe, very few previous 
studies deal directly with plasma phenomena. The effect of a plasma in cosmology normally 
has been studied with respect to the origin of the magnetic field. For example, the study of 
Harrison |]J elaborates a model of the origin of the magnetic field due to turbulence in the 
primordial plasma. There have been studies that analyze the effect of a magnetic field 
on primordial nucleosynthesis. Other studies, like that of Halcomb et al. ||, deduced the 
dispersion relation of waves, taking into account the expansion of the universe. There still 
lacks a general study of plasma phenomena related to cosmology. 

A plasma in thermal equilibrium, sustains fluctuations of the magnetic field (even for 
a non-magnetized plasma). This study is concerned with the study of the magnetic field 
spectrum in the primordial plasma. 

The electromagnetic fluctuations in a plasma has been made in numerous works, includ- 
ing those of Dawson ||, Rostoker et al. ||, Sitenko et al. ||, and Akhiezer J7|. Most of 
the results are compiled in the books of Sitenko and Akhiezer et al. |8|||. Little attention 
has been given to the question of how the magnetic field spectrum looks in a plasma. A 
naive answer to this question might be that it is a blackbody spectrum with a cut-off at the 
plasma frequency, knowing that photons only propagate in a plasma for uj > u p , where u p 
is the plasma frequency. This is not true however, due to the magnetic fluctuations of the 
plasma. 

Cable and Tajima and Tajima et al. [|lOHl2| performed a broad study of the magnetic 
field fluctuations in a plasma. They based their analyses on the Fluctuation-Dissipation The- 
orem. They were concerned, in particular, with the low-frequency spectrum of fluctuations, 
because, as they pointed out, no expression exists. 

The Fluctuation- Dissipation Theorem p,^] predicts the intensity of electromagnetic fluc- 
tuations. The intensity of such fluctuations is highly dependent on how the plasma is 
described, in particular, on the dissipation mechanisms used. 



3 



Cable and Tajima and Tajima et al. [|T0HT2[] studied the magnetic field fluctuations, 



for several cases. Two of their descriptions concern the primordial plasma which we are 
interested in, which is an isotropic, non-magnetized and non-degenerate plasma: a) a cold, 
gaseous plasma and b) a warm, gaseous plasma described by kinetic theory. 

In their study, Cable and Tajima (hereafter CT) in case (a) used the cold plasma de- 
scription with a constant collision frequency. In case (b) they analyzed the spectrum of 
fluctuations only for low frequencies, with the warm plasma description for phase velocity 
uj/k less or equal to the thermal velocity of the electrons, v e and the ions, Vi in a collisionless 
description. 

For the cold plasma description the spectrum that they obtain has a large zero-frequency 
peak. As the frequency is increased, the spectrum first drops below the blackbody spectrum 
in vacuum, then becomes the blackbody spectrum at high frequencies. In case (b), for the 
warm plasma description, the analyses was made only for the low frequency regime and 
they argued that the zero-frequency peak is present as well. They argue that the energy 
contained in the peak is approximately equal to the energy lost by the plasma cut-off effect. 

In order to obtain a correct magnetic field spectrum, it is necessary to describe the 
plasma in the most complete way as possible, taking into account thermal and collisional 
effects in a unified description. 

In this study we extend the pioneering work of CT, presenting a model that includes, in 
a unified description, collisional and thermal effects. Our model is based on kinetic theory 
incorporating thermal effects for all frequencies and wave numbers (not only for uj/k < 
v e , Vi). In order to describe the collisions that exist in the plasma, we used a model collision 
term. This collision term describes binary collisions, as used in the work by CT. In this way, 
we extend the previous model describing thermal and collisonal effects for all frequencies 
and wave numbers. Their description, the cold plasma and the warm plasma description 
in the collisionless case, are special cases of this model. (A pure collisionless treatment is 
unreal, such (b) extended to all frequencies, since if there were no collisions, then 

only Cherenkov emission could produce fluctuations, and there are no particles traveling 



fast enough to emit light waves.) 

However, for a fully ionized plasma as is our case, a treatment that takes into account 
collisions in a more complete way is necessary. Our model, an extention of the CT model, 
describes the basic features of a kinetic description. 

We present in Section |TJ the general expressions of the magnetic field fluctuations based 
on the Fluctuation- Dissipation Theorem. We review the cold plasma description in Section 
[I A| and the warm plasma description in the collisionless case, in Section [II B| . In Section [II C 



we present a general discussion and criticism of the assumptions made by CT. In Section 



UJ we present our model. Finally, in Section [TV], we discuss the results and present our 
conclusions. 

II. MAGNETIC FIELD FLUCTUATIONS 

The spectrum of fluctuations of the electric field in a plasma, given by the Fluctuation- 
Dissipation Theorem (for the deduction of the Fluctuation- Dissipation Theorem from the 
general relation of fluctuations in a plasma see ||[|) is, 

where 

Aij(u, k) = ^ (^f - 8i\ + £ij (u, k) , (2.2) 
where Eij(u, k) is the dielectric tensor of the plasma. For an isotropic plasma, 

A« = ^+(*«-f)(--^), (2.3) 

where El and Et are, respectively, the longitudinal and transverse dielectric permittivities 
of the plasma. In this case ||, 

h I hkj Im Et / kik~\ Im Et 

<^>w = ^SB^tI j^TJjr + («» " ii) ~ Z,-f \ ■ {2A) 



E T - 
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Using the fact that Bkw = (c/w)k x E^, we have the expression for the magnetic field 
fluctuations, 

, „ ft ( r kiki\ ( kc\ 2 Im et ,„ _ N 



£T- (77 



So, 



(2.6) 



An intuitive way to understand the above expression, is that the Fluctuation- Dissipation 
Theorem takes into account the emission and absorption processes in a plasma, and knowing 
that in equilibrium they are equal, the fluctuation level is evaluated. As Tajima et al. (1992) 
point out Hl2| , "an individual mode decays by a certain dissipation, giving up energy to 



particles or other modes, while particles (or other modes) excite new modes and repeat the 
process and the amount of fluctuations is related to the dissipation" . 

In short, to determine the fluctuations of the magnetic field in a plasma in equilib- 
rium or quasi-equilibrium, it is sufficient to know the transverse dielectric permittivity of 
the plasma, in particular, the dissipation mechanisms present for each frequency and wave 
number (Im Et)- This depends on the treatment used to describe the plasma. 

Another important feature of the magnetic field spectral distribution, can be seen from 
Eq. (|2.6|) . The equation 

e T (uj,k)-(^J = (2.7) 

determines the transverse eigenfrequencies of the plasma. Therefore, the magnetic field 
spectral distribution has steep maximae at the frequencies that correspond to the transverse 
plasma eigenfrequencies. In the transparency region (Im Et *C Re et), the magnetic spec- 
trum has 5-function-like maximae near the eigenfrequencies (i.e., the frequency spectrum of 
the fluctuations contains only the transverse eigenfrequencies in the plasma, Knowing 
that the photons have the dispersion relation in the plasma u 2 = lo 2 p + k 2 c 2 , we note that 
for frequencies uo ^> uo p the eigenfrequencies manifest themselves and the magnetic field 
spectrum behaves like a blackbody spectrum. 



A. Cold Plasma Description 



The cold plasma description does not take into account the thermal movement of the 
electrons, describing them by fluid equations. In such a description, collisionless damping 
is not included. To include collisions in this case is straightforward, being only necessary to 
add a new term oc rjv in the fluid equations, where r\ is the collision frequency and v the 
velocity of the particles. We follow here the CT study. 

CT used a multifluid model of the plasma (neglecting the v x B forces due to the small- 
ness of the velocities and the electromagnetic fields), 

ma ~dT = 6 ~ r]amaVa ' ( ' 

where a is a particle species label and rj a is the collision frequency of species a. From 
the above equation (performing a Fourier transformation and re-arranging the terms), the 
dielectric tensor can be obtained: 

CT studied the case of an electron-positron plasma. This is the plasma that dominated the 
universe in the beginning of primordial nucleosynthesis, at T ~ 1 MeV. In the case of an 
electron-positron plasma, u 2 e+ = u 2 e ^ and i] e + = rj e - = r\. Eq. ( |2.9| ) then becomes 

uj 2 

eij(u,k) = 5ij — — -5ij , (2.10) 

uo{uo + ir]) 

where u p 2 = u) pe + 2 + u) pe - 2 . r/ e was taken as the Coulomb collision frequency, r\ e = 2.91 x 
10" 6 n e lnAT-*l 2 (eV) s' 1 , where n e is the electron (positron) density. (For the case of an 
electron-proton plasma, also treated by CT, rj p = 4.78 x 10~ 18 n e lnAT~ 3 ^ 2 (eV) s^ 1 ). That 
is, this collision frequency describes the binary collisions in a plasma. 
For an isotropic plasma, the transverse dielectric permittivity is, 

*r(tt,k) = l- U l. . (2.11) 
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Substituting Eq. (|2.11|) in Eq. ( |2.6|) , CT obtained the magnetic field spectrum in the 
cold plasma description. If relativistic temperatures effects are included, the substitution 
Up — ► Up/ v /7 is made. 

The magnetic field spectrum (B 2 )^ is found by integrating (B 2 ) ku over wave numbers 
k (and dividing by (27r) 3 ). The (B 2 )^ diverges for high wave numbers. CT deal with this 
problem, breaking the integration on wave number into two intervals. One interval runs 
from | k |= to | k |= k cut . The other interval runs from | k |= k cut to | k |= oo. In the first 
interval, they keep the collision frequency r\ finite and in the second interval they let r\ — ► 
and drop the low-frequency part of the spectrum. The final expression obtained was 



87T ~ TT 2 e (fcWT)a/ - 1 '\cj Jo (J 2 + ^ 2 )x 4 + ... 



2tj P= / dx 



2 ,2,3/2 



h(u>' — u. 
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^ ) — J/ ^ (^) e^-o;^] , (2.12) 



2tt fe( Ba W 2 V _ij 

where 9 is the Heaviside step function, J = u/ujp e , lu p ' = u p /u pe , rj' = r]/uo pe and x = kc/uo pe . 



The first term extends up to the frequency uJhev = y k 2 ut c + uo 2 . The second term is the high- 
frequency and high-wave number expression (i.e., the spectrum for frequencies u > uJhev)- 

The justification given to this break-up procedure was that 77 should vanish smoothly 
as k — > 00 and as long as the results "do not critically depend on the manner in which 
77 approaches zero" the abrupt cut-off should be acceptable as a crude model. No strong 
justification was given to the choice of k cu t- They chose x cu t = k cut c/uj pe = 1 mainly because 
this is the value that makes the frequency spectrum be smooth at the joining of the low- 
frequency and the blackbody spectrum plO. 



B. Warm Plasma Description 

The warm plasma description describes the plasma based on the kinetic theory that 
takes into account the thermal distribution of the particles. In this description collisionless 
damping, like Landau damping, appears. The kinetic theory is based on the BBGKY hi- 
erarchy equations that describe a system of many particles. These equations are solved by 
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expanding the distribution function of the many particles in terms of the plasma parameter 
g = l/n\\), where A^> is the Debye length and n is the particle density. The description 
that is used usually is the collisionless description, based on the Vlasov equation that does 
not take into account collisions (and neglecting v x B forces): 

\ot m J 

where V = <9/<9x and V v = d/dv. The Vlasov equation in first order (in g) takes into 
account collisions, where the term on the right hand side of the equation is now the collision 
term: 

(! + T . V + £ B .V v )/(*,v,«)=(g) o . (2.14, 

For the collisionless description (assuming an isotropic plasma), the transverse dielectric 
permittivity is obtained from Eq. Q2.13 ). Assuming that the particles have a Maxwellian 



velocity distribution, the transverse dielectric permittivity is [|l3|-|i~6|, 

vClVv 2 [ , ( UJ \ . / 'ty\ 1 / 2 uj ( uj 2 



^•^-vfrn^rw ^{-^j}' (2 - i5) 

where a is the label for each species of the plasma and v a is the thermal velocity for each 



species, v a = ^T/m a . z = u/\^2kv a and 4>(z) = 2ze~ z Jq e x dx. 

In order to investigate thermal effects, CT used the collisionless description (i.e., the 
Vlasov equation in the regime where cu/k is less than the thermal speed of the plasma 
constituents). They studied a hydrogen plasma, thus the region investigated was u/k < v e , Vi 



where v e = yT/m and Vi = yT/M. (The extension to an electron-positron plasma, as in 
case (a), is straightforward: M — > m so v j — >• v e +). They used, therefore, Eq. (|2.15|) 



in the limit z < 1 and expanded the plasma dispersion function 4>(z). Substituting the 
approximate expression of ^^(cUjk) (expanded in the limit of z 1) in Eq. fl2.6|) , they 
obtained (B 2 ) kLU /8n. 

CT examined only the low frequency behavior of (B 2 )^, noting that it diverges for high 
wave numbers and a cut-off procedure is necessary. They assumed the same upper limit as 

in case (a), x U pp er x cu t = 1. 



C. Discussion and Criticisms 



For the cold plasma description, used in case (a), CT obtained a large zero- frequency 
peak in the magnetic field fluctuation spectrum (arguing that the warm plasma description 
exhibits it as well). 

The total spectrum is obtained with the link of the low frequency term to the high 
frequency term (Eq. ( |2.12 )). The behavior of the spectrum can be seen in Figure la. We 
study an electron-positron plasma at T = 7 x 10 9 K and n e = 4.6 x 10 30 cm -3 . This is 
the plasma in the beginning of primordial nucleosynthesis. The magnetic field spectrum 
S(u) = (£> 2 ) w /87r is divided by a normalization So = u^ksT '/c 3 . The dashed curve is the 
first term of Eq. (12.12 ), the low frequency spectrum. The dash-dot-dash curve is the second 
term of Eq. ( |2.12|) , the high frequency spectrum (obtained with rj — > 0). We note the link 
point used between the two curves. The blackbody spectrum in vacuum is also plotted 
(solid curve). It can be seen that the general behavior is that after the peak, with increasing 
frequency, the spectrum drops below the blackbody spectrum. At high frequencies, it merges 
into the blackbody spectrum. In Figure lb we plotted only the cold plasma spectrum for 
low and high frequencies. In Figure lc is a electron-proton plasma at T = 10 9 K and 
n e = 5.4 x 10 26 . In the epoch of primordial nucleosynthesis, at lower temperatures, the 
electrons and positrons annihilate and the plasma is reduced to a plasma of protons and 
electrons. In Figure Id we used the same plasma as in Figure lc, but plotted only the cold 
plasma spectrum for low and high frequencies. 

CT argued that the results do not critically depend on the upper limit. We show below 
that this is not true. As CT noted, the divergence occurs due to the subtle interaction 
between matter and radiation in small scales. CT used a classical fluid equation with a con- 
stant collision frequency (Coulomb collision frequency). This collision frequency describes 
the binary interactions in the plasma. They chose the cut-off x cu t = 1, basically because this 
is the value that makes the frequency spectrum smooth at the joining of the low-frequency 
and the blackbody spectrum. 
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In section VII of CT, they gave a quantum-mechanical justification of not extending k 
beyond k cut . They argue that for (hk) 2 /2m 3> k^T the plasma has a negligible effect on the 
electromagnetic spectrum. Let us call this k, ku m ((hku m ) 2 / 2m = ksT). 

When treating Coulomb collisions, a cut-off has to be taken, since at small distances 
the energy of the Coulomb interactions of the particles exceeds their kinetic energy which 
violates the applicability of the condition of the perturbation expansion (in the plasma 
parameter g <C 1). This occurs approximately for distances r m j n ~ e 2 /T or, more exactly, 
for the distance of closest approach between a test particle and an electron in a plasma, 
kmax = 1/ 1 f'min — Mmv 2 / {m + M) | eq \, where M, v and q are respectively, the mass, 
velocity and charge of the test particle [15]. 

Comparing the value of k cut , ku m and k max we see that k cut (x cut = k cut c/uj pe = 1) is much 
smaller than the others. For example, for the cases that CT used: 1) T = 10 10 K, n e = 
4.8 x 10 30 cm- 3 ; 2) T = 10 6 K, n e = 6.5 x 10 9 cm" 3 ; and 3) T = 10 4 K, n e = 6.5 x 10 3 cm" 3 
we found: 1) x Um = k Hm c/u pe = 11.5 and x max = k max c/u pe = 2444.4; 2) x Um = 3.1 x 10 9 
and x max = 6.6 x 10 9 ; and 3) xu m = 3.1 x 10 10 and x max = 6.6 x 10 11 . 

As x cut is quite arbitrary, we decided to vary the upper limit in the first term in Eq. 
( |2.12j ), until x upper = Xmax to observe the behaviour of the curves. This is shown in Figure 
2 where again we plot S(u)/Sq vs u/uj pe . 

We plotted the low frequency spectrum term not only until u hev = \Jk 2 cut c 2 + uj p 2 , as 
CT did, but extended it to higher frequencies. The dash-dot-dash curve has the upper limit 
Xupper = x cut as CT used. The dotted curve has the upper limit x upper = 2 x cut ] the dashed 
curve has the upper limit x upper = 5 x cut ; and the dash-two dot-dash curve has the upper 
limit x upper = x max . We also plotted the blackbody spectrum in vacuum (solid curve). It can 
be seen that by extending the upper limit to higher and higher values we obtain more and 
more of the blackbody spectrum, before the spectrum drops. This is easy to understand 
since, as we noted before, the magnetic field spectrum, obtained from the Fluctuation- 
Dissipation Theorem, contains the transverse eigenfrequencies of the plasma, the photons. 
When we extend the value of the upper limit, we permit higher eigenfrequencies to manifest 
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themselves. This occurs up to the frequency when, due to the upper limit chosen, no more 
photons can manifest themselves. 

Another feature that appears is that by increasing the upper limit the valley that appears, 
due to the plasma cut-off effect, becomes smaller and smaller. Eventually, for x upper = x max , 
we obtain a large peak for uj ~ but without any valley. In fact, for frequencies u < 2u pe 
the curve is above the blackbody spectrum in vacuum. This can be seen in Figure 3a, where 
the dashed curve is the low frequency term (in Eq. (|2.12|) ) extended to high frequencies (not 
only for u> < uJhev) with x upper = x max compared to the blackbody spectrum in vacuum (solid 
curve). In Figure 3b we plotted only the cold plasma curve, extended to high frequencies, 
showing that with x upper = x max the full blackbody spectrum is reproduced. 

CT argued, based on the behavior of the cold plasma spectrum with x upper = x cut , that 
the energy under the ui ~ peak is approximately equal to the energy stolen from the 
blackbody spectrum due to the plasma cut-off. In their words, this happens because, the 
"plasma squeezes the fluctuation energy of modes with frequency less than uj p into modes 
with frequencies very close to zero" . 

First of all, clearly, this does not happen if the upper limit is high enough in order to 
reproduce the full blackbody spectrum at high frequencies (with x upper = x max ) because the 
entire curve is above the blackbody spectrum. Second, we noted previously that a blackbody 
(photon) spectrum with a cut-off uj < uj p is expected taking into account that the photons 
only propagate for u > uo p in a plasma, where uo p is the plasma frequency. However, when 
we speak about blackbody, we speak about modes which propagate, in our case, photons. 
In vacuum, the photons have a dispersion relation u = kc and they are present in the entire 
frequency spectrum. In a plasma, when the dispersion relation is u 2 = uj 2 + k 2 c 2 , they only 
appear for to > uj p . 

Dawson 0] deduced the radiation spectrum in a plasma by a Gedanken experiment. 
A slab of plasma at a temperature T is put between two blackbodies at temperature T. 
Between the plasma and the blackbodies are vacuum regions. Radiation is emitted by the 
blackbodies and enters the plasma. In equilibrium, the plasma radiates the same amount of 
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radiation that it absorbs. Dawson deduced the density of radiation in a plasma as: 




(2.16) 



where U p is the radiation density in the plasma and U v the blackbody radiation density in 
the vacuum. It can be seen that a new factor appears, due to the presence of the plasma. 
This kind of Gedanken experiment tells us, however, only about modes that propagate. 
Nothing is told about the modes that do not propagate which appear due to correlations 
in the plasma. They are present at low frequencies, uj < u> p , but they also are present at 
uj > uj p , contributing, besides the photons, to the magnetic field spectrum. Our results show 
that this happens only for very high frequencies (u ^> uj p ) when the photons dominate the 
magnetic field spectrum. The magnetic field fluctuation spectrum has little to do with the 
photon blackbody spectrum in vacuum for u < uj p . In principle, the magnetic field spectrum 
can be greater, or less, than the photon blackbody spectrum in vacuum for u < uj p . The only 
manner to obtain the magnetic field spectrum is analyzing the magnetic field fluctuations 
from, for example, the Fluctuation-Dissipation Theorem. Using the cold plasma description, 
for example, with the cut-off x upper ~ x max , we find that the magnetic spectrum has more 
energy than the blackbody photon spectrum in vacuum. 



III. OUR MODEL 

In order to extend the work of CT and have a more complete description of the plasma, 
we desire a model that includes thermal effects as well as collisional effects. For this, we 
need a kinetic description that takes into account collisions. We used the Vlasov equation 
in first order. This equation gives collisional corrections to the Vlasov equation in zero 
order (collisionless case). It retains terms of order g and neglects terms of higher order. 
As the plasma parameter g is much less than unity (for example, as pointed out by |12 



in the epoch of t = 10 2 — 10° s in the primordial universe, g = 10 3 ), this is a good 
approximation, as terms of higher order are much smaller. As noted before, the term on the 
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right hand side of Equation (2.14), (df/dt)c, is the collision term. Obtaining (df /dt)c is a 
matter of great difficulty and different forms are required for various types of collisions, such 
as electron-electron, electron-neutral molecule, etc. Most of the expressions for (df/dt)c 
involve integral functionals of the distribution function /. The basic difficulty in taking into 
account the effect of collisions lies in the complexity of the solution of the kinetic equation 
when the correct collision integral is used. The problem can be simplified if a model collision 
term is used, that is, an approximate expression. 

The Boltzmann collision term takes into account all the possible binary collisions which 
the particle under observation might suffer. It is applied to weakly ionized plasma, when 
the scattering of charged particles by neutrals is predominant. In a fully ionized plasma (as 
in our case), the collisions are not predominantly binary and (df/dt)c for a test particle 
does not derive mainly from the possibility that other particles approach very closely and 
abruptly deflect it. The cumulative effect of more distant particles is more important. The 
charged particles simultaneously interact with all the other particles in the Debye sphere, 
which is a large number. The Fokker-Planck collision term is appropriate for a fully ionized 
plasma. It is based in the notion that a large-angle deflection of a particle by collisions is 
produced more rapidly by a succession of small-angle scattering with distant particles. It 
takes into account the effect of microscopic fields produced by all the other particles in the 
plasma. A test particle then is subject to simultaneously "grazing" collisions and its progress 
in velocity space becomes a random walk. 

We used the BGK collision term as a rough guide to the inclusion of collisions in the 
plasma. BGK is a model equation of the Boltzmann collision term. (For a derivation see 



Clemnow et al. |16] and Alexandrov et al. [0.) As CT, we are then only treating binary 



collisions. (Concerning this collision term, see for example the work of Sitenko and Gurin 
H, who studied the effect of an effective binary collision frequency on the fluctuations in a 
plasma). A more complete treatment using the Fokker-Planck collision term is necessary. 
However, due to its complexity, as it involves integral functionals of the distribution function 
/, the kinetic equations are extremely difficult to solve. We used as an effective collision 
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frequency, the Coulomb collision frequency as CT. 

The BGK collision term conserves the number of particles. The universe at high temper- 
atures, at the beginning of the nucleosynthesis (T ~ 1 MeV) is an electron-positron plasma. 
As it cools down (T < 0.5 MeV) the electron and positrons start to annihilate and finally, 
at low temperatures, the plasma is reduced to a plasma of protons and electrons. Let us see, 
if this is a good approximation in the early universe. The rate of annihilation is given by 
r ~ nva, where n is the density, v the average velocity and a the average cross section. The 
average cross section is given by, a = ^[ (3 ~^ )A -2 + /3 2 ], where (3 2 = 1 — A = ln[j^\, 
a = and s = (2E) 2 with m the electron mass and E the electron energy [18|]. Comparing 
T with the collision frequency rj (the Coulomb collision frequency as used by CT), we see 
that r < 77. For example, for T ~ 0.8 MeV, V = 4 x 10 16 s' 1 and r] = 4 x 10 17 s -1 . 
Therefore, the characteristic time of the kinetic processes is shorter than for annihilation 
and the approximation of the conservation of the number of particles at high temperatures is 
valid. (As in our case, when the collision frequency is much less than the plasma frequency, 
the collision time is the dominant time for the kinetic processes.) At low temperatures, an 
electron-proton plasma has to be considered. 

We note that because the electron-positron plasma strictly does not conserve the particle 
number due to annihilation and creation, there may not be a great advantage of using the 
BGK collision term, although the above paragraph indicates that the annhilation frequency 
is very much smaller than the collision frequency. 

The plasma, in the epoch of primordial nucleosynthesis, is almost collisionless. (For 
example, for T = 10 10 K, r]/uj pe = 10~ 3 .) Therefore, a possible procedure that one might 
think of is to expand the dielectric permittivity in terms of rj/co. However, we are interested 
in the spectrum for all frequencies, even to ~ 0, and no matter how small 77 is, we require 
frequencies with u < 77. Thus an expansion in rj/u cannot be made. 

Therefore, we perform an analyses of the Vlasov equation in first order with the BGK 
collision term, without making any approximation, 
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'df 

dt 



= ~V(f ~ /max) , (3.1) 

C 



where rj is the Coulomb collision frequency (considered constant). f max is given by 
/max(x, v, t) = iV(x, t) f (v) / N where the number density is N = N + Ni and / = f + fi, 
fo being the unperturbed Maxwellian distribution. Substituting this collision term in the 
equation of Vlasov in first order, Eq. ( ^.14[ ), performing a Fourier transformation and re- 
arranging the terms, we arrive at 



/i(v) ' 



i(u — irj — k ■ v) 



e„dfo N ! / 
Ei 7T + ^^T^° 



(3.2) 



To eliminate N\ we integrate over velocity space JIB 



N, = J fi(v)dv. (3.3) 

We have the current j = e J vfi(v)dv. Knowing that = (JijEj, then for an isotropic plasma 
the transverse permittivity is easily obtained (generalizing for several species): 

where a is the label for each specie of the plasma, v a is the thermal velocity for each specie 
and Z(z) is the Fried & Conte function [Il9l, 



y7T J-oo r — Z 

If relativistic temperatures effects are included, the substitution u p — > u p / ^7 is made. 

It can be seen that no approximation needs to be made on io/\/2kv a and 77/ \/2kv a . (Re- 
writing Z(z) in terms of the error function, Z(z) = i^/7rexp(— ^ 2 )[1 + erf(i^)] [|19j, it can be 
easily solved numerically without having to take any asymptotic limit). 

Another thing that we need to be careful of is related to the dielectric permittivities in 
the region of u/k > c, where c is the velocity of light. For u/k > c, the term connected with 
Cherenkov emission in the imaginary part of the dielectric permittivities has to vanish. This 
is because there are no particles with velocities greater than the speed of light to produce 
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Cherenkov emission. (The condition for Cherenkov emission of a wave u(k) requires particles 
with velocities satisfying the condition v — ou/k = v p , where v p is the phase velocity of the 
wave.) For our case, a plasma in thermal equilibrium, we assume as CT in section V of 
their article [|11|], that the distribution function is a Maxwellian one. The calculation of 
the permittivity tensor (and the dielectric permittivities) involves integral functionals of the 
distribution function. As pointed out by Melrose ||14|| , the non- vanishing of the imaginary 
parts, in a collisionless treatment (that is Cherenkov emission) is due entirely to the effects of 
unphysical particles with v > c in the Maxwellian distribution. We deal with this problem by 
requiring that the Cherenkov emission in the imaginary part is zero in the regime u/k > c. 
The term of the imaginary part, connected with the collisional damping, is not set to zero. 

As photons have phase velocities greater than the speed of light (u/k > c), the Cherenkov 
emission cannot produce the photons. In a pure collisionless treatment, the imaginary part 
has to be set to zero in the regime of u/k > c. In this case, no photons are produced and the 
treatment is unrealistic. In our model, which includes collisions and thermal effects, only 
the term connected with the Cherenkov emission has to be set to zero. In this treatment 
we ensure that the photons are produced only by collisions, that is, by bremsstrahlung 
(free-free) emission. 

Therefore, by forcing the Cherenkov emission in the imaginary part to be zero in the 
regime of u/k > c, we ensure that no spurious effects contaminate the result. We do this 
not only in the regime of photons (u > u pe ), but for all frequencies. That is, Cherenkov 
emission of magnetic fluctuations occur only in physical regimes, i.e., for u/k < c. It is 
left for a future study a fully relativistic treatment (needed for the primordial universe at 
high temperatures). There, the imaginary parts of the dielectric permittivities are zero for 

> c, where v$ = u/k, when 7 = (1 — v 2 /c 2 ) and p^ = "y^mv^ are imaginary. 

In principle, the fact that the term connected with the Cherenkov emission, for u/k > c is 
forced to be zero, could be troublesome. This procedure, however, did not change appreciably 
the results. For example, the change of the intensity of the magnetic field fluctuations 
A(B 2 ) uj (after integrating on wave number) for T = 7 x 10 9 K, n e = 4.6 x 10 30 cm -3 is: 
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10~ 7 , 10~ 3 and IO -5 for uo/uj p = 0.1, 1.0 and 100.0, respectively. The change is the difference 
between (ImET) c ^ and (ImET) c = 0, where (Ime^Q is the term in the imaginary part 
of the transverse dielectric permittivity connected with Cherenkov emission. 

Before we go on, it is interesting to comment and emphasize that both the cold plasma 
description and the warm plasma collisionless description are particular solutions of this 
model. Let us consider, as CT did, an electron-positron plasma. For | z | 3> 1, where 
z = (u + irj)/\/2kv e , we expect to obtain the cold plasma approximation. In the limit of 
| z | 2 3> 1, Z(z) = — - — ^3 + ••• — iy/nze~ z2 . The last term is due to Cherenkov emission, 
but as we noted before, in this limit this is a spurious solution and we need to set it to 



zero. Taking only the first term, and substituting in Eq. (|3.4j ), we obtain the cold plasma 
dielectric permittivity, as expected. 

The warm plasma collisionless description is obtained by setting rj ^ 0. In this limit, if 
we write z = x + iy, iy = iO, 4>(x + iO) = </>(x) — iy/iixe~ x2 (Z(z) = —1/z 0), Eq. ( |3.4|) is 
equal to the warm plasma collisionless dielectric permittivity, as expected. 

We substitute the dielectric permittivity Eq. fl3.4Q in Eq. ( ^.6| ) obtaining the magnetic 
field spectrum {B 2 ) ku) . Here, a divergence at high wave numbers also occurs. Let us discuss 
this in detail. Our model uses a kinetic theory description with a model collision term that 
describes the binary collisions in the plasma. In our case, a cut-off has to be taken, since for 
very small distances the energy of the Coulomb interactions of the particles exceeds their 
kinetic energy which violates the applicability of the condition of the perturbation expansion 
(in the plasma parameter g 1). This occurs approximately for distances r min ~ e 2 /T, or 
more exactly, the distance of closest approach between a test particle and an electron in a 
plasma, k max = l/r min = Mmv 2 /{m + M) \ eq |, where M, v and q are respectively, the 



mass, velocity and charge of the test particle ]K 

The cut-off procedure can only be removed, treating properly the effects of distant en- 
counters. This can be done with the Fokker-Planck collision term. Thompson and Hubbard, 
and Hubbard in several works pp|-F2^|, analyzed the Fokker-Planck equation and its coef- 



ficients. The diffusion and friction coefficients that appear in the Fokker-Planck equation, 



take into account correlation effects between distant particles in the plasma. When higher 
order terms in the Fokker-Planck equation are calculated and summed, a term resembling 
the Boltzmann collision term is obtained. In their treatment, they showed that the cut-off 
procedure is unnecessary. However, due to the complexity of the solution of the kinetic 
equation with the Fokker-Planck collision term, we used the BGK collision term. With this 
model collision term, a cut-off is necessary and we chose x max consistent with this model 
collision term. A more exact treatment, however, is needed. 

In Figure 4a we plot the magnetic field spectrum S(u) = (5 2 ) w /87r (divided by the 
normalization So = w^ksT/c 3 ) vs uj/uj pe for an electron-positron plasma at T = 7 x 10 9 K 
and n e = 4.6 x 10 30 cm -3 . The dotted curve is our model and we compare it with the 
blackbody spectrum in vacuum (the solid curve). In Figure 4b, we extend the curves to high 
frequencies, showing the behaviour of the blackbody at high frequencies. Figures 4c and 4d 
are for an electron-proton plasma, with T = 10 9 K and n e = 5.4 x 10 26 . 

As we commented before (section |11 C| ), the results are very dependent on the cut-off 
chosen. Using x max as the cut-off, we obtain results that differ from the CT results: First, 
the peak intensity found by them for frequencies u ~ decreases. This is due to the kinetic 
plasma effects that smear out the peak. However, it is interesting to see that qualitative 
agreement between the work of CT and ours exist with respect to the zero frequency peak; 
Second, we obtain the blackbody naturally for high frequencies; and Third, the magnetic 
field spectrum has more energy than the blackbody spectrum for frequencies u < 4 u pe . 



IV. CONCLUSIONS AND DISCUSSION 

The magnetic field spectrum can be deduced from the fluctuations of the magnetic field 
described by the Fluctuation- Dissipation Theorem and it is highly dependent on the way 
the plasma is described. We discussed the cold plasma description and the warm plasma 
description in the collisionless case studied by CT. We showed that x cu t is much smaller than 
x max and xi im , where xu m = k Um c/u) pe used in Sec. VII of CT (arguing that for k > ku m the 
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plasma has a negligble effect). x max is the cut-off used in treating binary collisions, (xm^)" 1 
being the distance of closest approach between a test particle and an electron in a plasma 
(divided by c/u pe ). 

We also showed that the Fluctuation- Dissipation Theorem contains the eigenfrequencies 
of the plasma (in the transverse case), the photons. Using the cold plasma description 
with the upper limit x upper = x max , we obtain the blackbody spectrum at high frequencies 
naturally, without the necessity of a link procedure used by CT. For this case, the valley 
disappears and the curve is above the blackbody spectrum in vacuum. 

The calculations were made for two types of primordial plasmas: The electron-positron 
plasma at the beginning of the Big Bang nucleosynthesis; and the electron-proton plasma 
at lower temperatures. 

The manner to obtain the entire magnetic field spectrum is analyzing the magnetic field 
fluctuations. This is the only way to obtain information, not only about modes that propa- 
gate (i.e., photons), but also modes that do not propagate. The modes that do not propagate 
appear not only at low frequencies but also at high frequencies due to the correlations in the 
plasma. Only at very high frequencies does the photon contribution dominate the magnetic 
field spectrum. The argument used by CT, that the energy under the peak is almost equal 
to the energy "stolen" by the plasma cut-off effect of the blackbody in vacuum, is incor- 
rect. There is no reason why we have to have the same energy as the blackbody spectrum 
in vacuum for photons for u < u pe , since the photons have a different dispersion relation 
than the fluctuations in the plasma. In fact, using a upper limit x upper = x max in the cold 
plasma description, for example, the spectrum obtained is above the blackbody spectrum in 
vacuum. 

The reason why the collective modes of the plasma can have more energy for uo < uo p 
than the photons in vacuum, can be understood as follows. Photons are massless bosons 
with the dispersion relation oj 2 = k 2 c 2 . For the energy interval, < uj < u p , the wave 
number interval is k = to k equal to u p /c. A relatively small amount of phase space is 
involved. For the collective motions of the plasma, in general, we have a larger amount of 
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phase space. For example, for plasmons with energy uo ~ uj p , the amount of phase space 
extends to a maximum k of hp = u> p /v t , which is greater than uj p /c for the photons. In 
general, for a given frequency for uj < u p , the greater phase available to the collective modes 
of the plasma (than that of the photons) implies more energy, or a higher spectrum. 

We presented a model that incorporates, in the same description, thermal and collisional 
effects. We used the Vlasov equation with the BGK collision term. This collision term 
describes the binary collisions in the plasma. A model that takes into account collisions 
in a more complete way is necessary. For a fully ionized plasma it is necessary to use the 
Fokker-Planck collision term that takes into account the effect of the microscopic fields. Due 
to the complexity of the solution of the kinetic equation with such a collision term, we used 
the BGK collision term. This model, an extention of the CT model, describes the basic 
features of a kinetic description. 

As we noted before, the results are very dependent on the cut-off chosen. Using x max as 
the cut-off, consistent with the collision term used, we obtain results that differ from the CT 
results. The final magnetic spectrum of a non-magnetized plasma in thermal equilibrium 
has the following characteristics: a) The peak intensity found by CT for frequencies u ~ 
decreases; b) The blackbody is obtained naturally for high frequencies; and c) The magnetic 
spectrum has more energy than the blackbody photon spectrum in vacuum, in particular, 
for frequencies, u> < 4cj pe , where uj pe is the electron plasma frequency. 
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FIGURES 

FIG. 1. The magnetic field spectrum ln[S(u)/So] vs oo/u> pe (where S(uj) = {B 2 )^/^ and 
So = Lo pe 2 kBT/c 3 is the normalization) for the cold plasma description with x upper = x cu t ~ 1 for: 
(a) The electron-positron plasma at T = 7 x 10 9 K and n e = 4.6 x 10 30 cm" 3 (the dashed curve is 
the cold plasma spectrum for low frequencies, the dash-dot-dash curve is the cold plasma spectrum 
for high frequencies (the link point between the two curves is indicated), and the solid curve is 
the blackbody spectrum in vacuum); (b) The same as case (a), where the cold plasma spectrum is 
plotted extended to high frequencies; (c) The same as case (a) but for an electron-proton plasma 
at T = 10 9 K and n e = 5.4 x 10 26 ; and (d) The same as case (b) but for an electron-proton plasma 
at T = W 9 K and n e = 5.4 x 10 26 . 

FIG. 2. The magnetic field spectrum ln[S(uj)/So] vs oj/uj pe for the cold plasma description 
for various upper limits (where S(u) = {B 2 ) w /&-k and Sq = uj pe kBT/c 3 is the normalization). 
The plasma is an electron-positron plasma at T = 7 x 10 9 K and n e = 4.6 x 10 30 cm" 3 . The 
dash-dot-dash curve is the spectrum with x uvper = x cu t', the dotted curve is the spectrum with 
Xupper = 2 x cut ; the dashed curve is the spectrum with x upper = 5 x cu t; and the dash- two dot-dash 
curve is the spectrum with x upper = x max . The solid curve is the blackbody spectrum in vacuum. 

FIG. 3. The magnetic field spectrum ln[S(u)/So] vs io/to pe (where S(u>) = (B 2 ) uJ /87r and 
So = uj pe kBT/c 3 is the normalization) for the cold plasma description with x upper = x max for: (a) 
The electron-positron plasma at T = 7 x 10 9 K and n e = 4.6 x 10 30 cm~ 3 (the dashed curve is 
the cold plasma spectrum and the solid curve is the blackbody spectrum in vacuum); and (b) The 
same as case (a), extended to high frequencies showing the blackbody behaviour. 
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FIG. 4. The magnetic field spectrum ln[S(bj)/So] vs uo/uJ pe (where S{oj) = {B 2 ) w /Stt and 
Sq = ujp e kBT/c 3 is the normalization) for our model for: (a) The electron-positron plasma at 
T = 7 x 10 9 K and n e = 4.6 x 10 30 cm~ 3 (the dotted curve is the spectrum of our model and 
the solid curve is the blackbody spectrum in vacuum); (b) The same as case (a), extended to 
high frequencies; (c) The same as case (a) but for an electron-proton plasma at T = W 9 K and 
n e = 5.4 x 10 26 ; and (d) The same as case (c) extended to high frequencies. 
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